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^ . Abstract 

H ' 

' We consider the production of neutrinos in active galactic nuclei 



(AGN) through hadronic cascades. The initial, high energy nucleons 
are accelerated in a source above the accretion disk around the central 
black hole. From the source, the particles diffuse back to the disk and 
initiate hadronic cascades. The observable output from the cascade 
are electromagnetic radiation and neutrinos. We use the observed 
diffuse background X-ray luminosity, which presumably results from 
this process, to predict the diffuse neutrino flux close to existing limits 
from the Frejus experiment. The resulting neutrino spectrum is E~ 2 
down to the GeV region. We discuss modifications of this scenario 
which reduce the predicted neutrino flux. 
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1 Introduction 



The active nuclei of galaxies (AGN), ranging in luminosity from Seyfert galax- 
ies to quasars, are the most powerful individual sources of radiation in the 
Universe. To explain the power emitted by such objects, one generally as- 
sumes the existence of a central engine in which the gravitational energy 
of matter falling into a supermassive black hole gets converted into radia- 
tion. Even though so far only electromagnetic radiation has been observed, 
it is generally assumed that other particles are accelerated and emitted as 
well to explain the tight relationship between the non-thermal and thermal 
components in the UV and X-ray spectra [I], |3|, |, Hj. Of special inter- 
est are neutrinos, since they can travel cosmological distances without losing 
the information on the direction they originated from. Large underwater 
detectors || |7j or detectors in the antarctic ice cap || are used as neutrino 
telescopes. Recent calculations have shown that the flux of neutrinos origi- 



nating in an AGN could be detected by such experiments |9|, [II], [TTJ, [12], [L3|]. 
Data from proton decay experiments [14| and airshower arrays [15] is already 
sensitive enough to constrain such models significantly. 



2 The AGN model 

In the following we are interested in the production of neutrinos in radio- 
quiet AGN. In the "standard" model for AGN one assumes that the central 
black hole is surrounded by an accretion disk of infalling matter. Besides 
that, one expects to find bipolar outflow of gas and plasma perpendicular to 
the disk (jet). Jets can be seen in different objects with disk accretion and in 



many AGN |Tj| [L7|]. One expects shocks in the plasma of the jets |TI| which 
could accelerate protons through first order Fermi acceleration to energies 
up to O(10 9 GeV) with a powerlaw spectrum of E~ 2 [18|. The observation 



of 7-rays with the same kind of spectrum indicates that indeed there 



is shock-acceleration outside the core of the AGN, so that the photons can 
escape. 

In |[20|| , Niemeyer showed that the far infra-red (FIR) emission of AGN 



can be explained by assuming that the accelerated protons diffuse back from 
the jet to the accretion disk and heat dust clouds beyond the outer region 
of the accretion disk. In the inner part of the disk, protons hitting the disk 
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will initiate hadronic cascades through interactions with the gas in the disk. 
This could feed into an electromagnetic cascade and thereby generate the 
observed X-ray and gamma emission [T3| . 

To describe the accretion disk, we use the model by Shakura and Sun- 
For our purposes, we will concentrate on the inner region, which 
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yaev 

is the radiation dominated part of the disk (neglecting the thin, innermost 
ring around the black hole). 



3 Hadronic cascades 

Due to the observed short timescales of the variability of the X-ray and 
UV components of the AGN spectrum, this part of the electro-magnetic 
radiation must predominantly originate in a small, central region around the 
central black hole [fT0| , |9|] with a typical mass of 1O 8 M for luminous AGN. 
We expect neutrino production to take place in the same region, so we need 



to know the incoming particle flux in this region. From |2D we know that 
the FIR spectra are best fitted by assuming that the protons originate from 
a point source in the jet at zq < O(10 18 cm) m 3 ■ 10 4 i?s above the disk. We 
consider two possibilities how protons from this or a similar source can travel 
back to the central region of the disk and initiate hadronic cascades in the 
accreting gas. 

I) The region interesting for the production of neutrinos is much closer 
to the centre of the disk; it extends to about 200-Rs ~ 10 16 cm. In the disk 



model this corresponds to the innermost region of the disk. An isotropic 
proton source in the jet, similar to the kind of source mentioned above, 
located close to the origin, i.e., at z Q <C 200-Rs, could then deposit at least 
half its luminosity into the disk. This happens simply for geometrical reasons, 
since the inner disk occupies about half of the horizon of the proton source. 
In this approximation, there is no dependence on the exact distribution of 
the proton source along the jet or on details of the proton transport. 

II) It is also possible that the proton source at zo = O(10 18 cm), needed 
for the explanation of the FIR radiation, emits a strong downward proton 
flux. Such an anisotropy in the proton emission arises naturally, since the 
interaction of the protons in the source with upward UV to X-ray photons 
from the accretion disk preferably produces downward neutrons through the 
reaction p'-f —> A —> nii + , if the photon scattering in the ambient medium 
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is negligible. The upcoming photon density is sufficiently low that r„ 7 < 1, 
i.e., the neutrons will reach the disk without further reactions with photons. 
Furthermore, above 50 TeV, the neutrons will reach the disk without decaying 
back into protons. The incident nucleons will then induce a hadronic cascade 
in the inner accretion disk. In this scenario, we get beaming for high energy 
neutrinos. At lower energies, where protons are part of the cascade, the 
magnetic field in the disk randomises the orientation of the cascade. Also, 
the spectrum of the nucleons emitted from the source is flatter than E~ 2 at 
high energies which leads to a flatter high energy neutrino spectrum as well. 

In the following, we will concentrate on scenario I with a pure E~ 2 spec- 
trum for the proton source. In the inner region of the accretion disk, the 
column density of the disk is 

UW \3R S J cm 2 1 ; 

If we take M/M e dd ~ 0.1 and a ~ 0.1, we see that the column density is 
high enough compared to the mean free path for pp-collisions of O(50 g / cm 2 ) 
for a hadronic cascade to develop. The magnetic field confines the protons 
to the disk, which leads to an even further increase of the effective column 
density seen by the protons. Even though neutrons which are produced in 
the cascade are not confined by the magnetic field, the amount of gas present 
is sufficient to prevent them from escaping. 

The hadronic cascade in the disk is much simpler than cascades in the 
earth's atmosphere [^, |23], ^4j| , since the density of the accretion disk is 
so low that all unstable particles decay rather than interact. Therefore the 
cascade consists of a nucleonic part which feeds into the mesonic and electro- 
magnetic channels; no pion-nucleon reactions occur. Furthermore, the pion 
channel will always be the dominant channel for the production of neutrinos. 
This is different from the production of cosmic ray neutrinos in the atmo- 
sphere, since there reactions of pions with nuclei remove the pions from the 
parent population for neutrino production. This leads to a steeper neutrino 
spectrum and to the increasing importance of kaons and charmed mesons for 
the neutrino production in the atmosphere while heavy mesons are of little 
importance for AGN neutrino production. Under the assumptions of scaling 
hadronic cross-sections, the power law of the incident protons is kept and 
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Figure 1: The structure of the hadronic cascade: pp and pn interactions drive 
the cascade and feed the production of pions, leptons and hard photons. We 
do not distinguish particles and anti-particles in this figure. So the leptonic 
final channels need to be doubled. 



inherited by all the secondaries, i.e., the flux of particle x is 

n x (E) = C x h p (E), n p (E) = C E^~\ (2) 

where C is the normalisation of the proton spectrum. The specific power in 
protons is then 

-i(max) 



^ t0t) = Xd) E'n p (E')dE', 



= < 



C'(7-1)" 1 (^S) _7+1 , for 7 7^1, ® 

/ „(max) \ 

Cm f^rj, for 7=1- 



For the secondaries, the energy range is shifted down by a factor e x ; the 
resulting power is 



^ = /"(!, E'n p (E')dE', 



— C f-"/ +1 R( tot ) 



(4) 



The cascade equations can be solved approximately to determine C x 
and e x . For scaling hadronic cross-sections, the experimental data is sum- 
marised by the spectrum weighted moments 

Zab = I x^F^xl) dx L , (5) 
Jo 
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Species 




c x = 4 tot V^ tot) 


7 


0.070 


0.29 


e 


0.056 


0.16 


Vp from 7i 


0.040 


0.13 


from \i 


0.056 


0.16 




0.048 


0.14 



Table 1: Relative normalization of the spectrum of secondaries, shift of par- 
ticle energies, and fraction of the incident proton luminosity energy carried 
by the stable secondaries for 7=1. The total energy in one species is split 
evenly between particles and anti-particles. 



where xl = E b jE a is the fraction of the energy of the primary particle trans- 
ferred to the secondary and F ah {xj)) = E b dn b (xi) /dE b is the dimensionless, 
inclusive cross-section for the production of particles of type b. For 7 > 1, the 
weighting factor xj^ 1 reduces the importance of the small-x^ region, where 
the accuracy of experimental data is low. Unfortunately, in the case relevant 
here, there is no such supression of the small-x^ region. The percentage of 
energy not accounted for in the last column of table [I] is a rough measure 
for the uncertainty resulting from this. We see that ~ 49% of the power 
is emitted as neutrinos where ~ 33% are muon type neutrinos. The rest is 
electromagnetic where twice as much luminosity is emitted in photons as in 
electrons. Muon neutrinos and anti-neutrinos, which are the species which 
experiments are most sensitive to, carry about 2/3 of the electromagnetic 
power. 



3.1 Electromagnetic radiation 



The electromagnetic output of the hadronic cascade is reprocessed in the 
inner disk. Pair cascades, inverse Compton scattering, and reflection on cold 
material change the shape of the initial _E -2 -spectrum and provide a steep 

§, m, ||, |(J. This produces 



511 keV 25 



turnover around 

the observed X-ray and gamma emission of the AGN which we will calculate 
in detail elsewhere. 
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4 Resulting neutrino spectrum 



The neutrino spectrum from a single source mirrors the -E _2 -spectrum of 
the protons, only that it is shifted down by a factor of t v ~ 0.05. The 
upper cutoff of the neutrino spectrum depends on the details of the shock 
acceleration process in the jet, since that determines the maximum energy 



reached by the protons [13 



The Neutrino luminosity: To get an estimate of the neutrino luminosity 
of a source, we use its total emission in X-rays and 7-rays as a reference. We 
assume that the power emitted at X-ray and 7-ray energies is the electro- 
magnetic power output of the hadronic cascade, reprocessed in the inner disk 
to yield the observed spectrum (see section |3.1| ). The relation we get is then 

E^ = 0.95 4*$. (6) 

Using an observation for 4+7 > equation (£|) allows us to determine C v . There 
is actually a logarithmic dependence on the energy range of the neutrinos; 
for definiteness we take the range 10~ 2 GeV . . . 10 6 GeV. 



4.1 X-ray background 

To be able to use the relation between electromagnetic output and neutrino 
production to predict the diffuse neutrino background, we have to estimate 
how much of the observed, diffuse X-ray and 7-ray background results from 
hadronic cascades in AGN. There is growing evidence from the ROSAT all 
sky survey that an appreciable fraction of the background is due to active 



galactic nuclei, possibly a dominant proportion [ pip . 

The active nucleus in a quasar is often surrounded by a zone of very 
active star formation. ROSAT X-ray observations of the Seyfert galaxy 
NGC1068 demonstrate that even the 2-10 keV X-ray emission might 
be dominated by starburst activity near to the nucleus. On the other hand, 
X-ray observation of the starburst galaxy M82 show that the X-ray spec- 
trum is rather hard and is likely to contain an inverse Compton component, 
i.e., a powerlaw component. Hence the starburst surrounding the active 
nucleus is likely to contribute a hard powerlaw X-ray emission component. 
What fraction of the total emission from a radioweak quasar might then 
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arise from this hard X-ray emission component? In starburst galaxies the ra- 
dio (5 GHz), mid-infrared (60 /mi), and X-ray (2keV) emission is correlated 
(see || for a comparison of starburst galaxies with quasars), and thus the 
weakness of the radioemission attributable to the starburst speaks against a 
dominant contribution from the starburst region, as does the often observed 
variability of the X-rays emitted by radioweak quasars. On the other hand, 
this argument is not sufficiently general to eliminate a strong contribution 
by a starburst in all quasars. 

Thus we cannot decide conclusively to what extent a starburst surround- 
ing the active nucleus in many radioweak quasars contributes hard X-ray 
emission similar to that of the nuclear emission. As a matter of fact, ex- 
tensive work on compton reflection models for AGN X-ray emission cannot 



explain all features of the observed X-ray background spectrum |34j. To be 
conservative, we allow for a factor of three maximum between the total hard 
X-ray emission from an AGN galaxy (starburst and active nucleus together) 
and the contribution strictly from the hadronic cascade. 

The observations of the hard X-ray background at those photon energies 
minimally influenced by reprocessing ||, i.e., at energies where the original 
E~ 2 powerlaw is still visible, give a possible range for the unreprocessed 
energy total of 1.0- 10 5 to 1.4-10 5 eV cm -2 s^ 1 sr _1 . Using the abovementioned 
conservative estimate, we arrive at a lower limit of 3 • 10 4 eV cm -2 s _1 sr~ x for 
the contribution of hadronic cascades to the X-ray and 7-ray background. 

4.2 Neutrino background 

Since the cosmological redshift is the same for neutrinos and photons, we can 
scale the background neutrino luminosity using the fraction of the X-ray and 
7-ray background derived above. Using equations (|J) and @, we get 

N(E U ) = 1.6 • 1(T 12 (^r^J 2 cm" 2 s^ 1 sr- 1 GeV- 1 (7) 

as a conservative limit for the sum of all neutrino species. About 2/3 are 
muon neutrinos, the remaining 1/3 are electron neutrinos. This prediction 
is a factor of 2.5 lower than the experimental limit set by the Frejus experi- 



ment [14], For each family, the number splits evenly into neutrinos and 
anti-neutrinos. Again, we have a weak dependence of the scale factor on 
the range of the neutrino energies. This expression is valid up to energies 
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of ~ 0.05-E^, max ) = 10 7±1 GeV; the spectrum turns over at this energy due to 
the lack of parent particles to produce neutrinos. Actually, we expect the 
background spectrum to have a less sharp cutoff at high energies than a single 
source, since the maximum neutrino energy varies between different AGN. 
Therefore, the cutoff in the superposition of all spectra will be smoothed out. 

Such a spectrum is similar to the results of the Monte-Carlo calculation 
in jy]; but the flux of muon neutrinos predicted here is about one order of 
magnitude lower. Compared to the prediction for the neutrino background 
made in ||, we expect the E~ 2 powerlaw for the neutrino spectrum down to 
the GeV range without the flattening seen in || below 10 6 GeV — a difference 
which should be visible in the experiment. 



5 Discussion 

The main point in our model is the production of neutrinos as the result 
of a hadronic cascade. Compared to the p7-channel, which has a threshold 
of E p « 8 • 10 6 GeV for photoproduction on UV photons, all protons above 
a few hundred MeV contribute in pp-interactions. As a consequence, we 
expect neutrino emission even from sources with a low cutoff in the primary 
proton spectrum. This way, we include contributions from a larger class of 
sources both to the neutrino background and to the diffuse X-ray background. 



Similar to the result of other authors || [TjJ , our model displays the tendency 
to produce neutrinos strikingly close to existing limits [14, 35, 15 1. 

Modifications of this pp-model can lead to a reduced neutrino flux while 
reprocessing of the electromagnetic component ensures an unchanged X-ray 
emission. The most drastic modification — reducing the maximum proton 
energy and correspondingly the maximum neutrino energy — can ultimately 
decrease the observable extragalactic neutrino flux by moving the cutoff be- 
low the cross-over with the steeper spectrum of atmospheric neutrinos. More 
realistically, steepening of the AGN proton distribution at a sufficiently low 
break energy is an alternative to reduce the predicted event rate in a neu- 
trino detector. In both cases, neutrino production via p7-reactions becomes 
ineffective. 

On the other hand, a proton spectrum flatter than E~ 2 leads to the 
dominance of jry-reactions at high energies. This, again keeping the X-ray 
background unchanged, enhances the neutrino flux far above the TeV range 
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while reducing the flux below. Details depend on the modelling of the target 
photon field (i.e., contributions from the disk, corona, jet, etc.). 

Another important point in our model is that the acceleration of the 
protons takes place above the disk. This way, the acceleration takes place in 
an environment more suitable than an accretion disk. 
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